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Abstract 

We calculate the total gravitational energy and the gravitational energy den- 
sity of the de Sitter space using the definition of localized gravitational energy that 
naturally arises in the framework of the teleparallel equivalent of general relativity. 
We find that the gravitational energy can only be defined within the cosmological 
horizon and is largely concentrated in regions far from the center of spherical sym- 
metry, i.e., in the vicinity of the maximal spacelike radial coordinate R = \J\- The 
smaller the cosmological constant, the farther the concentration of energy. This re- 
sult complies with the phenomenological features of the de Sitter solution, namely, 
the existence of a radial acceleration directed away from the center of symmetry 
experienced by a test particle in the de Sitter space. Einstein already contemplated 
the de Sitter solution as a world with a surface distribution of matter, a picture 
which is in agreement with the present analysis. 
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I. Introduction 

The difficulty in defining gravitational energy density within the framework of the Hilbert- 
Einstein Lagrangian formulation has led to the belief that the gravitational energy cannot 
be localized. It is widely assumed that an expression for the localized energy density of 
the gravitational field does not exist. However it is well known that the total energy of 
a given asymptotically flat spacetime can be calculated by means of pseudotensor meth- 
ods, which make use of coordinate dependent expressions. A different approach to the 
construction of an energy expression for the gravitational fleld is based on the idea of 
quasilocal energy. The quasilocal deflnition of energy, momentum and angular momen- 
tum associates these quantities to an arbitrary spacelike two surface S in an arbitrary 
spacetime manifold. The inexistence of an unequivocal deflniton of gravitational energy 
still remains an actual problem, important in its own right. Furthermore such deflniton 
may play a major role in the thermodynamics of self-gravitating systems. This problem 
has been recently addressed in ref . |^ , where a comprehensive bibliography on quasi-local 
energy is presented. Although all attempts so far have led to interesting mathematical 
developments, they did not allow the achievement of a deflnite solution, either because of 
conceptual or mathematical difficulties. 

Recently the problem of localization of energy in general relativity has been recon- 
sidered in the framework of the teleparallel equivalent of general relativity (TEGR)p|. 
The Lagrangian formulation of the TEGR is established by means of the tetrad fleld e° ^ 
and the spin affine connection ujnab, which are taken to be completely independent fleld 
variables, even at the level of fleld equations. This formulation has been investigated 
in the past in the context of Poincare gauge theories^, ^. However, as we will explain 
ahead, this is not an alternative theory of gravity. This is just an alternative formulation 
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of general relativity, in which the curvature tensor constructed out of Ufj_ab vanishes, but 
the torsion tensor is non-vanishing. The physical content of the theory is dictated by 
Einstein's equations. As we will show, in this alternative geometrical formulation the 
gravitational energy density can be naturally defined. 

The expression for the localized energy density of the gravitational field has arisen in 
the context of the Hamiltonian formulation of the TEGRj5[. It has been demonstrated 
that under a suitable gauge fixing of cj^afe, already at the Lagrangian level, the Hamiltonian 
formulation of the TEGR is well defined^. The resulting constraints are first class 
constraints. In fact the Hamiltonian formulation looks very much similar to the the 
usual ADM formulation]^. However there are crucial differences. The integral form of 
the Hamiltonian constraint equation C = in the TEGR can be written in the form 
C = H — Eadm = 0, when we restrict considerations to asymptotically flat spacetimesP]. 
The quantity e{x) which appears in the expression of C and which under intergration 
yields Eadm is recognized as the gravitational energy density. We have applied the 
expression of e{x) to the calculation of the energy inside a surface of constant radius, 
both for the Schwarzschild]^] and for the Kerr metric]^, and the results are remarkably 
the same as those obtained by means of the quasilocal energy definition proposed by Brown 
and York|l[]. Moreover, the calculational scheme is rather simple, as we will see shortly, 
and is exempt of some complications inherent to the latter. Therefore for asymptotically 
fiat spacetimes the gravitational energy density has a definite and unambiguous expression 
within the framework of the TEGR. 

We recall that the gravitational energy can also be calculated by means of the surface 
term that appears in the expression of the gravitational Hamiltonian[^ |^. However, 
such surface term yields only the total gravitational energy, as the integration has to be 
necessarily carried out over the whole three dimensional spacelike hypersurface, in which 
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case the lapse function N{x) goes over into its asymptotic value ^ 1 at spatial infinity. 

The action integrals for spacetimes with different topologies require surface terms 
that are specific to each topology. Thus the corresponding Hamiltonian also acquires a 



surface term that is determined by the topological boundary conditions [|T21- However the 
Hamiltonian constraint for a spacetime foliated by spacelike hypersurfaces always has the 
same basic structure, irrespective of boundary conditions {additional terms such as the 
cosmological constant may appear in the Hamiltonian constraint, as we will see ahead in 
eq.(lO) ). 

Therefore the question immediately arises as to whether the Hamiltonian constraint 
equation in the TEGR can always be writtem as C = H — E = 0, in which case e{x) 
would be the gravitational energy density for any curved spacetime. 

One of the simplest deviations from asymptotically flat geometries are spacetimes 
with conical defects. We have applied our expression of gravitational energy density to 
the calculation of the energy per unit length of defects of topological nature, which include 



disclinations, i.e., cosmic strings, and dislocations |[TT|1. The result is quite encouraging. 
We arrive at precisely the same well known expression for the energy per unit length of 
a cosmic string (not even multiplicative factors have to be adjusted). Moreover the total 
energy of a dislocation is zero, a result which is in close analogy with the statements of 
the theory of elasticity, which asserts that disclinations and dislocations are defects which 
require high energy and low energy, respectively. 

In this paper we consider the de Sitter space, which is topologically of the type x R. 
We restrict the considerations to the static region within the cosmological horizon (i.e., 
the region for which —goo > ) and calculate both the total energy and the distribution 
of energy along the radial direction. Again the result is rather remarkable. We will show 
that the cosmological constant induces a distribution of gravitational energy in such a 



way that the energy is largely concentrated in the peripheral region, i.e., in the vicinity 



of the maximal spacelike radial coordinate R = ^ j^. As we will show in section III, 
this picture is in total agreement with the phenomenological features of the de Sitter 
solution, and is as well in agreement with Einstein's belief, according to which the de 
Sitter's solution represents a spacetime with a surface distribution of matter [|l2l. This 



fact strongly supports the validity of our expression for the gravitational energy density 
and also represents a clear indication that the Hamiltonian constraint equation in the 
TEGR can be unambiguously interpreted as an energy equation of the type H — E = 0. 

Notation: spacetime indices and local Lorentz indices a,b,... run from to 3. 

In the 3+1 decomposition latin indices from the middle of the alphabet indicate space 
indices according to /i = 0,z, a = (0),(i). The tetrad field e"^ and the spin con- 
nection ujfj^ab yield the usual definitions of the torsion and curvature tensors: i?" b/iu = 
df_iU>^"' b + ^fi"" c^u'^ b — T°- iiu = d^e"" y + oj^"" u — ■■■■ The flat spacetime metric is 
fixed by ?^(o){o) = -1- 



II. The Lagrangian and Hamiltonian formulations of the TEGR 

In the TEGR the tetrad field e°' ^ and the spin connection uo^ab are independent field 
variables, not related by any of the field equations. The spin connection is enforced to 
satisfy the condition of zero curvature. The Lagrangian density in empty spacetime is 
given by[|, I 

L(e,a;,A) = -fce( Jt-^^^T,,, + ^T'^^^T,,, - T^TJ + eA«''^'^i?,fe^,(u;) . (1) 
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where k = ys^' ^ gravitational constant; e = det{e°- Lagrange 

multipliers and is the trace of the torsion tensor defined by = 

The equivalence of the TEGR with Einstein's general relativity is guaranteed by the 
identity 

eR{e,u) = eR{e) + e{-T^''T,,, + ^T^^^T,,, - T^TJ - 2d^{eT^^) , (2) 

which is obtained by just substituting the arbitrary spin connection u^ab = °W/ia6(e) + 
Kfj^ab in the scalar curvature tensor R{e,uj) in the left hand side of (2); °u;^afe(e) is the 
Levi-Civita connection and K^ab = ^^b^iTx^y + Tj,a^ — T^^x) is the contorsion tensor. 
The vanishing of i?" bpLui'^)-, which is one of the field equations derived from (1), implies the 
equivalence of the scalar curvature -R(e), constructed out of only, and the quadratic 
combination of the torsion tensor. It also ensures that the field equation arising from the 
variation of L with respect to e"" ^ is strictly equivalent to Einstein's equations in tetrad 
form. Let = denote the field equation satisfied by Ca^. It can be shown by explicit 
calculations that 

1 r„ 1 

(we refer the reader to ref. for additional details). 

For asymptoticaly flat spacetimes the total divergence in (2) does not contribute to 
the action integral. Therefore the latter does not require additional surface terms, as it is 
already invariant under coordinate transformations that preserve the asymptotic structure 
of the field quantities^]. It is well known that for compact geometries a surface term has 
to be included in the action, in order to make the variations of the field variables well 
defined. This surface term is constructed out of the trace of the extrinsic curvature on 
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the boundary. However we will no longer worry about surface terms in the Lagrangian or 
in the Hamiltonian, as we will be interested only in the constraint structure of the theory. 

The Hamiltonian formulation of the TEGR can be successfully implemented if we fix 
the gauge uoab = from the outset, since in this case the constraints constitute a first 
class set|^. The condition Woafe = is achieved by breaking the local Lorentz symmetry 
of (1). We still make use of the residual time independent gauge symmetry to fix the 
usual time gauge condition e^k) ° = e(o)i = 0. Because of Woafe = 0, H does not depend 
on p''"-^^ the momentum canonically conjugated to cokab- Therefore arbitrary variations of 
L = pq — H with respect to P^"'^ yields cOkab = 0. Thus in view of uoab = 0, u>kab drops 
out from our considerations. The above gauge fixing can be understood as the fixation of 
a global reference frame. 

Under the above gauge fixing the canonical action integral obtained from (1) becomes]^ 

Atl = J d'x{H(^>e(,), - H} , (3) 

H = NC + N'Ci + S^^„H™" . (4) 

In expression (4) above we are omitting surface terms. and A^* are the lapse and shift 
functions, H™" = cq) "^H^-^^" and S^n = — are Lagrange multipliers. The constraints 
are defined by 

C = dj{2keT^) - keE'^'^n^j - ^(n^'n,, - ^H^) , (5) 

Cfc = -eo>a,H(^')^ - H(^')*r(,),, , (6) 

with e = det{e(j)k) and = g'^e^^'^^T(^j)ik, T(^j)ik = die(^j)k - dkC^y. We remark that (3) 
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and (4) are invariant under global SO (3) and general coordinate transformations (in eqs. 
(1) and (2) e is the determinant of the spacetime tetrad field; from eq.(3) on e stands for 
the determinant of the triads restricted to the three dimensional spacelike hypersurface) . 

If we assume the asymptotic behaviour e(j)fc ~ 'qjk + ^hjk{^) for r — oo, which is 
appropriate for an asymptotically flat spacetime, then in view of the relation 

/ d'^xdjieT^) = / dSk{dihik-dkh,,) = Eadm (7) 

where the surface integral is evaluated for r ^ oo, we note that the integral form of the 
Hamiltonian constraint C = may be rewritten as 

j d'x^keT.^'^n,^ + 4^(n*'n,. - ^n2)| = Eadm • (8) 

The integration is over the whole three dimensional space. Given that dj{eT^) is a scalar 
density, from (7) and (8) we define the gravitational energy density enclosed by a volume 
V of the space as[3 



Eg = 



^— J^d'xd,{eT^) . (9) 



It must be noted that this expression is also invariant under global SO (3) transformations. 

One is immediately led to ask whether the Hamiltonian constraint for topologically 
different spacetimes can also be written as eq.(8). In the next section we will consider 
the de Sitter space. Before addressing the latter, let us recall here some applications 
of Eg. We have calculated the gravitational energy inside a surface of constant radius 
To both for the Schwarzschild]^] and for the Kerr solution^, using Boyer and Lindquist 



coordinates 1 13, 14]. These quantities have also been calculated by means of Brown and 



York's precedure, in refs.[0] and [|T^, respectively. The expressions found by using (9) 
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are in total agreement with those obtained via the method of ref.[0]. Moreover Eg can 
be calculated for any volume in the three-dimensional spacelike hypersurface, as least 



through numerical integration, whereas the evaluation of the energy in ref.[jT5[ can only 
be carried out in the limit of slow rotation of the black hole (the application of Brown 
and York's procedure to the Kerr solution with arbitrary parameters meets some technical 
difficulties, as discussed in ref.[0). 

Definition (9) has also been applied to a class of conical spacetime defects, in which 
disclinations (cosmic strings) and dislocations are considered altogether. For the space- 
time of a single cosmic string, i.e., for a pure disclination, we obtain precisely the well 



known value of energy per unit length of the string|Tl|. Furthermore the total gravita- 
tional energy for a pure dislocation vanishes. This is a very interesting result, because 
we know from the theory of elasticity that disclinations are defects that require a large 
ammount of energy to be formed, whereas dislocations require low energy (see sections 



6.3.2 and 6.5 of ref. [|T6| for a discussion as to why the energy demanded for the formation 
of a disclination in a crystal is very high). Therefore the above calculations of energy are 
in close agreement with the statements of the theory of elasticity (in this respect we recall 
that attempts were made long time ago which envisaged the spacetime as a continuum 
with microstructure (see[§|, section 1.2)). 



III. Gravitational Energy in the de Sitter Space 

We will consider now the theory defined by the Lagrangian density (1) supplemented by 
a term containing the cosmological constant A. Thus we add to (1) the quantity 2'^eA, 
where = Ne is the determinant of the spacetime tetrad field e^^. This additional 
term will contribute to the action integral (3) only as an extra term of the Hamiltonian 
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constraint. The new Hamiltonian constraint reads 



C = d^{2keT^) - keE'^'Tkij - ^(^'"'^ji - ^n^) - 2eA , (10) 

The most general spherically symmetric solution of the field equations with a positive 
cosmological constant is the Schwarzschild-de Sitter solution (throughout this section we 
will make G — 1): 



ds' ^ -(l- — - ^)dt' + (l-'^-^) 'dr' + r'de' + r'stn^edcp^ . (11) 
V r R^J V r R^J 

This metric represents the gravitational field of a particle of mass m located at the origin 
of a globally hyperbolic spacetime. The vacuum solution, obtained by setting m — 
in (11), is the de Sitter solution. R is the maximal spacelike radial coordinate for the 
(vacuum) de Sitter space and is given hy R = y^. 

Strictly speaking de Sitter spacetime corresponds to a four-dimensional surface in a 
fiat five-dimensional space with metric (-,-|-,-|-,-|-,-|-) described by 

-z^ + zf + zi + zl + zl + zl = J, A > . 

The coordinates {t,r,,9^(f)) in (11) cover only half of the space defined by the relation 
above. However we will be interested just in (11), as it suffices for our purposes. Moreover 
we will restrict the considerations to the physical region between the Schwarzschild (black 
hole) and the cosmological horizons. 

Expression (9) allows us to calculate the gravitational energy for any volume in space. 
We wish to obtain the energy contained within a surface of constant radius Tq. For this 
purpose we will calculate eT^ — eT^ for a spacetime whose spacelike section is described 
by the line element 
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(12) 



where a is a function of the coordinate r. The triads that correspond to (12) are given 

by 



/ a sin6cos(f) r cos9cos(f) —r sinOsincf) \ 

^{k)i — a sinOsincf) rcosOsincj) r sin6cos(j) . (13) 

\ a cos9 —r sin9 / 

{k) is the hne index, and i is the column index. 

The determinant e of (13) reads e = ar^sinO. After a lengthy but otherwise straight- 
forward calculation of 



(k)ij 



we arrive at 



eT^ = 2rsine(l-^^ 



Therefore for the Schwarzschild-de Sitter solution we have 



(14) 



The energy contained within a surface of constant radius is thus given by 



(15) 



where is a surface of constant radius Tq. 

Let us evaluate expression (16) for the range of values of To such that 



(16) 
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2m Tf? 

■ 

in which case we assume the cosmological constant to be very small. Expanding (16) and 
neglecting all powers of both ^ and ^ we arrive at 

Eg = ro + m = Eds + m . (17) 

Efis is the value of energy we would obtain in the absence of the mass m, and therefore 
it represents the background (vacuum) energy. Upon subtraction of this term we obtain 
the standard ADM value of energy for a particle of mass m. Of course in (17) we expect 
To to be much larger than m. 

The total gravitational energy contained within the cosmological horizon can be easily 
calculated, but for this purpose one has to find the roots of the equation 1 — ^ — =0. 
The result is not illuminating. If i? ^ m we find that = r^ax, where Vmax is 

slightly smaller than R and is a solution of the equation above. In what follows we will 
rather consider the vacuum de Sitter solution only, since in this case the analysis of the 
energy density is most easily carried out, and the main features are not altered by the 
introduction of a mass m at r = 0. 

Before proceeding we mention that the present analysis is different from that carried 



out by Abbott and Deser[|T^. These authors provide an expression for the energy of the 
gravitational field about the de Sitter background, i.e., they calculate the energy of a 
field configuration that deviates from the de Sitter metric and which vanish at infinity. In 
contrast, by means of expression (9) we can compute the energy of the whole gravitational 
field configuration, including the background. 

The total gravitational energy E^is contained in the physical region of the vacuum de 
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Sitter space is obtained from (16) by making m — and = R: 



EdS ^ R = ^j. (18) 

The total volume of the compact spacelike section equals 27r^i?^. Therefore the average 
energy density is given by 



Eas A ^^9) 



27rm^ 67r2 ' 

Let us obtain now the distribution of gravitational energy in the de Sitter space. In view 
of the spherical symmetry we will be interested in calculating the density of energy per 
unit radial distance s{r), which is obtained by integrating -^dr{eT^) in 9 and 0. Thus 
£{r) yields the gravitational energy contained between the spherical shells of radii r and 
r + dr. Considering m = in (15) we obtain upon integration in the angular variables 
and differentiation in r: 

2/32-1 

where we have set P'^ — The function e{r) has the following properties. In the range 
< /3 < 1 e(r) vanishes only for r — 0. Moreover for /3 = 1 it diverges: e{l3 = 1) ^ oo. 
It is straightforward to check that this is a monotonically increasing function, largely 
concentrated in the vicinity of /3 = 1: e(/3 = 0.1) = 0.015 ; e{P = 0.5) = 0.423 ; e{P = 
0.9) = 2,422. The total energy contained inside the surfaces of radii O.IR , 0.5R , 0.9R 
are given by Eg = 5.01 x IQ-'^R , Eg^O, 067R , Eg = 0.51R, respectively. 

Therefore almost half of the gravitational energy is located between /3 = 0.9 and 
P — 1. This result is in total agreement with the phenomenological features of the de 
Sitter solution, and can be verified in the following way. The goo component of (11) can 
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be written as 



goo = 1 + 2^ , 

where (j) is given by 

= Ar. 

r 6 

(f) is the potential in classical mechanics which would induce motion of a test particle 
approximately along the geodesies of (11). Therefore even in the absence of a mass m a 
test particle would be subject to a radial acceleration 

a = ^Ar , 

directed away from r = 0. 

The acceleration increases with the distance r, indicating that the gravitational field 
is more intense at points far from the origin. Therefore when m — the gravitational 
given by (11) acts on physical bodies as if there were a radially symmetric distribution of 
matter about the origin, beyond the cosmological horizon, just as m represents the mass 
of a black hole, concentrated inside the black hole horizon. 

This is precisely the picture we obtain from (20). By applying (9) to the de Sitter 
solution we find that the cosmological constant induces a spherically symmetric distribu- 
tion of gravitational energy, concentrated in regions distant from the origin, due to the 
gravitational field that acts on a test particle placed in the vacuum de Sitter space. Such 
a field can be thought as due to some matter distribution. 

One may think of (11) as representing the gravitational field of a spherical cavity inside 
some spherically symmetric distribution of matter (this idea is discussed, for instance, in 
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ref.|jT^). In this respect we recall that Einstein already conjectured that de Sitter solution 
would correspond to a world with a surface distribution of matter |T^. Such conjecture 



has found a natural explanation within the present geometrical framework, and shows 
that (9) yields a consistent expression for the gravitational energy in the de Sitter space. 

We will briefly discuss how our procedure applies to the anti-de Sitter solution. The 
latter is obtained by making the replacement — »■ in (11). The anti-de Sitter 
space is a non-compact manifold with constant negative curvature. The energy contained 
within a surface of constant radius r^, can be easily calculated and reads 



where we have ignored the mass term m. ranges from to oo. Therefore as — oo, 
we find that Eg — oo. This is an expected result, since the anti-de Sitter space is 
non-compact. The density of energy per unit radial distance e{r) in this case is given by 

1 + ^ 

We find that e{r) = only for r = 0. This point is also the only global maximum for e{r); 
for r — s> oo we clearly see that e{r) — oo. Thus e{r) is a non-positive monotonically 
decreasing function. 

IV. Discussion 

The definition of gravitational energy is a long-standing problem in the theory of general 
relativity. Numerous attempts have been made in the past for a solution. This problem 
still attracts considerable attention in the literature, and remains an important issue to be 
settled. Essentially all of these previous attempts are in one or another way unsatisfactory. 
In particular it is widely claimed that the gravitational energy cannot be localized. We 
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do not share this opinion. The mathematical structure of the TEGR shows that not 
only we do have a consistent and unambiguous definition of gravitational energy for 
asymptotically flat spacetimes, naturally built in the Hamiltonian formulation, but also 
that the gravitational energy is localized. The gravitational energy in the framework of 
the TEGR is given by expression (9). This expression has been successfully applied to a 
number of spacetimes, as we mentioned in section II, whose gravitational energy is already 
known. A justification for the extension of this definition to more general spacetimes is not 
straightforward. In the case of asymptotically flat spacetimes the Hamiltonian constraint 
equation can be written as C — H — Eadm — 0. We assume that this form of the 
constraint is a general feature of the theory, namely, that we can write the Hamiltonian 
constraint as C = 1/ — £^ for an arbitrary spacetime, since the constraint structure in 
general relativity is fixed and does not depend on any particular topology. 

In the above we considered the de Sitter solution and concluded that the cosmological 
constant induces a distribution of gravitational energy largely concentrated in the vicinity 
of the maximal spacelike radial distance R. This result is in total agreement with the fact 
that a test particle in the de Sitter space is subject to a radial acceleration directed away 
from the center of symmetry. Therefore the outcome of our analysis complies with the 
phenomenological behaviour of a test particle in the de Sitter space. To our knowledge 
this is the first time that such analysis has been provided. 

By inspecting equation (18) we see that if we make A — > the total energy Eg diverges. 
The vanishing of A in (11) ammounts to a change from a compact to a non-compact 
topology. Therefore such a change is not smooth, as it requires an infinite ammount of 
energy. This fact seems to indicate that, at the classical level, topology changing processes 
are forbidden. 



15 



Acknowledgements 

This work was supported in part by CNPQ. 



References 

[1] J. David Brown and J. W. York, Jr, Phys. Rev. D47, 1407 (1993) 
[2] J. W. Maluf, J. Math. Phys. 36, 4242 (1995) 

[3] F. W. Hehl, in Proc. of the 6th Course of the School of Cosmology and Cravitation 
on Spin, Torsion, Rotation and Supergravity, Erice, 1979, P. G. Bergmann and V. de 
Sabbata, eds., (Plenum, New York, 1980) 

[4] F. W. Hehl, J. D. McCrea, E. W. Mielke and Y. Ne'eman, Phys. Rep. 258, 1 (1995) 

[5] J. W. Maluf, J. Math. Phys. 35, 335 (1994) 

[6] R. Arnowitt, S. Deser and C. W. Misner, in Cravitation: an Introduction to Current 
Research, L. Witten, ed., (Wiley, New York, 1962) 

[7] J. W. Maluf and Andreas Kneip, Localization of Energy for a Kerr Black Hole, 
Universidade de Brasilia preprint (1995) ( |gr-qc/ 95040 11| ) 

[8] T. Regge and C. Teitelboim, Ann. Phys. 88, 286 (1974) 

[9] L. D. Faddeev, Sov. Phys. Usp. 25, 130 (1982) 

[10] S. W. Hawking and G. T. Horowitz, The Cravitational Hamiltonian, Action, Entropy 
and Surface Terms, preprint DAMTP/R 94-52, UCSBTH-94-37 ( fer-qc/95010T4D 



16 



[11] J. W Maluf and Andreas Kneip, Gravitational Energy of Conical Defects, Universi- 
dade de Brasilia preprint (1995) ( |gr-qc/9505008[ ) 

[12] W. Pauli, Theory of Relativity, (Dover, New York, 1958), section 62, page 182; A. 
Einstein, S. B. preuss. Akad. Wiss., 270 (1918); de Sitter, Proc. Acad. Sci. Amst. 20, 
1309 (1918) 

[13] R. P. Kerr, Phys. Rev. Lett. 11, 237 (1963) 

[14] R. H. Boyer and R. W. Lindquist, J. Math. Phys. 8, 265 (1967) 

[15] E. A. Martinez, Phys. Rev. D50, 4920 (1994) 

[16] N. Rivier, in Geometry in Condensed Matter Physics, Ed. J. F. Sadoc (World Scien- 
tific, Singapure, 1990) 

[17] L. F. Abbott and S. Deser, Nucl. Phys. B195, 76 (1982) 

[18] R. Adler, M. Bazin and M. Shiffer, Introduction to General Relativity, Second Edition 
(McGraw-Hill Kogakusha, Tokyo, 1965), page 399 



17 



